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Abstract. For a Levy process £ = (£,t)t>o drifting to — oo, we define the so-called exponential 
functional as follows 

f oo 

I f = / e^dt. 
Jo 

Under mild conditions on £, we show that the following factorization of exponential functionals 

h = 1 h- x ly 

holds, where, x stands for the product of independent random variables, H~ is the descending 
ladder height process of £ and Y is a spectrally positive Levy process with a negative mean 
constructed from its ascending ladder height process. As a by-product, we generate an integral 
or power series representation for the law of for a large class of Levy processes with two-sided 
jumps and also derive some new distributional properties. The proof of our main result relies 
on a fine Markovian study of a class of generalized Ornstein-Uhlenbeck processes which is of 
independent interest on its own. We use and refine an alternative approach of studying the 
stationary measure of a Markov process which avoids some technicalities and difficulties that 
appear in the classical method of employing the generator of the dual Markov process. 
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1. Introduction and main results 



We are interested in studying the law of the so-called exponential functional of Levy processes 
which is defined as follows 

/•oo 

L = / e^dt, 



where £ = (£t)t>o is a Levy process starting from and drifting to — oo. Recall that a Levy 
process £ is a process with stationary and independent increments and its law is characterized 
completely by its Levy-Khintchine exponent ^> which takes the following form 

r i a 2 f°° 

(1.1) logE e za =^( z ) = bz + —z 2 + (e zy - 1 - zyl {{yl<1} ) U(dy), for any z G iR, 

where a > 0, b G R and II is a Levy measure satisfying the condition j R (y 2 A l)H(dy) < oo. See 
PQ for more information on Levy processes. 

The exponential functional 1^ has attracted the interest of many researchers over the last 
two decades. This is mostly due to the prominent role played by the law of 1^ in the study of 
important processes, such as self-similar Markov processes, fragmentation and branching pro- 
cesses but also in various settings ranging from astrophysics, biology to financial and insurance 
mathematics, see the survey paper [6]. 

So far there are two main approaches which have been developed and used to derive infor- 
mation about the law of the exponential functional. The first one uses the fact that the Mellin 
transform of Lt is a solution to a functional equation, see (|4.ip below, and is due to Carmona 
et al. [TO] and has been extended by Maulik and Zwart [21]. It is important to note that (|4.ip 
is useful only under the additional assumption that £ possesses some finite, positive exponential 
moments since then it is defined on a strip in the complex plane. This equation can be solved 
for exponential functionals of negative of subordinators and spectrally positive Levy processes 
yielding some simple expressions for their positive and negative integer moments respectively, 
which, in both cases, determine the law. Recently, Kuznetsov and Pardo [16] have used some 
special instances of Levy processes, for which the solution of the functional equation can directly 
be guessed and verified from (|4.ip . to derive some information concerning the law of It. It is 
worth pointing out that, in general, it is not an easy exercise to invert the Mellin (or moments) 
transform of 1^ since a fine analysis of its asymptotic behavior is required. This Mellin transform 
approach relies on two difficult tasks: to find a solution of the functional equation and to provide 
a general criterion to ensure the uniqueness of its solution. For instance, this approach does not 
seem to successfully cope with the whole class of spectrally negative Levy processes. 

The second methodology, which has been developed recently by the second author in [23] 
and |26] , is based on the well-known relation between the law of Ig and the distribution of the 
absorption time of positive self-similar Markov processes which were introduced by Lamperti 
|18j in the context of limit theorems for Markov processes. Indeed, in |26| . it is shown that the 
law of l£ can be expressed as an invariant function of a transient Ornstein-Uhlenbeck companion 
process to the self-similar Markov process. Using some potential theoretical devices, a power 
series and a contour integral representation of the density is provided when £ is a possibly killed 
spectrally negative Levy process. 

In this paper, starting from a large class of Levy processes, we show that the law of Ig 
can be factorized into the product of independent exponential functionals associated with two 
companion Levy processes, namely the descending ladder height process of £ and a spectrally 
positive Levy process constructed from its ascending ladder height process. It is well-known 
that these two subordinators appear in the Wiener-Hopf factorization of Levy processes. The 
laws of these exponential functionals are uniquely determined either by their positive or negative 
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integer moments. Moreover, whenever the law of any of these can be expanded in series we can 
in general develop the law of in series. Thus, for example, the requirements put on the Levy 
measure of £ in [16] can be relaxed to conditions only on the positive jumps (the Levy measure 
on the positive half-line) of £ thus enlarging considerably the class of Levy processes £, for which 
we can obtain a series expansion of the law of 1^. 

Although our main result may have a formal explanation through the Wiener-Hopf factoriza- 
tion combined with the functional equation (|4.ip , the proof is rather complicated and involves a 
careful study of some generalized Ornstein-Uhlenbeck (for short GOU) processes, different from 
the ones mentioned above. For this purpose, we deepen a technique used by Carmona et al. [TU1 
Proposition 2.1] and further developed in [T7], which relates the law of 1^ to the stationary 
measure of a GOU process. More precisely, we show that the density function of 1^, say m^, is, 
under very mild conditions, the unique function satisfying the equation Cm^ = 0, where C is an 
"integrated infinitesimal" operator, which is strictly of an integral form. The latter allows for 
a smooth and effortless application of Mellin and Fourier transforms. We believe this method 
itself will attract some attention as it removes generic difficulties related to the study of the 
invariant measure via the dual Markov process such as the lack of smoothness properties for the 
density of the stationary measure and also application of transforms which usually requires the 
use of Fubini Theorem which is difficult to verify when dealing with non-local operators. 

Before stating our main result let us introduce some notation. First, since in our setting 
£ drifts to — oo, it is well-known that the ascending (resp. descending) ladder height process 
H + = (H + (t))t>o (resp. H~ = —H~ * = (—H~'*(t))t>o) is a killed (resp. proper) subordinator. 
Then, we write, for any z € iR, 

(1.2) <f>+(z) = logE [exp(zH+(l)j\ =5 + z+ [ (e 2 « - l)/x+(dy) - k + , 

J(0,oa) 

where <5+ > is the drift and k + > is the killing rate. Similarly, with <5_ > 0, we have 

(1.3) <f>-(z) =logE [exp(zH-(l))] = -6-z- f (1 - e _ ^)//_(dy) . 

We recall that the integrability condition Jq°°(1 A y)fJL±(dy) < oo holds. The Wiener-Hopf 
factorization then reads off as follows 

(1.4) *(z) = -c<j) + {z)(j)^{z) = -0 + (z)0_(z), for any z G iR, 

where we have used the convention that the local times have been normalized in a way that 
c = 1, see (5.3.1) in [13]. We avoid further discussion as we assume (|1.4p holds with c = 1. 

Definition 1.1. We denote by V the set of positive measures on M + which admit a non- 
increasing density. 

Before we formulate the main result of our paper we introduce the two main hypothesis: 

('Hi) Assume further that — oo < E [£i] and that one of the following conditions holds: 

E_)_: fi+ € V and there exists z+ > such that for all z with, 5R(z) E (0, z+), we have 

\$?(z)\ < oo. 
P+: 11+ € V. 
(H2) Assume that 

P±: /i+ € V, k + > and jt_ G V. 
Then the following result holds. 

Theorem 1.2. Assume that £ is a Levy process that drifts to —00 with characteristics of the 
ladder height processes as in (jl.2p and (jl.3p . Let either (Hi) or (H2) holds. Then, in both cases, 
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there exists a spectrally positive Levy process Y with a negative mean whose Laplace exponent 
ip+ takes the form 

POD 

(1.5) tp+{~ s ) = -s<M-s) = <M 2 + k +s + s 2 e~ sy n + (y, oo)dy, s > 0, 

J o 

and the following factorization holds 

(1.6) Ig = 1 H - X Iy 

where = stands for the identity in law and x for the product of independent random variables. 

Remark 1.3. We mention that the case when the mean is — oo together with other problems 
will be treated in a subsequent study as it demands techniques different from the spirit of this 
paper. 

The result in Theorem 11.21 can be looked at from another perspective. Let us have two 
subordinators with Levy measures n± such that fi + € V, k + > and /i_ € V . Then according 
to Vigon's theory of philanthropy, see [31], we can construct a process £ such that its ladder 
height processes have exponents as in (jl.2p and (11 .31) and hence £ satisfies the conditions of 
Theorem 11.21 Therefore we will be able to synthesize examples starting from the building 
blocks, i.e. the ladder height processes. We state this as a separate result. 

Corollary 1.4. Let \x± be the Levy measures of two subordinators and /x + 6 V, k + > and 
[i- € V . Then there exists a Levy process which drifts to —oo whose ascending and descending 
ladder height processes have the Laplace exponents respectively (|l,2p and (|1.3p . Then all the 
claims of Theorem \1.2\ hold and in particular we have the factorization (jl.6p . 

We postpone the proof of the Theorem to the Section [U In the next section, we provide some 
interesting consequences whose proofs will be given in Section [5j Finally, in Section [3l we state 
and prove several results concerning some generalized Ornstein-Uhlenbeck processes. They will 
be useful for our main proof and since they have an independent interest, we present them in a 
separate section. 

2. Some consequences of Theorem 11.21 

Theorem 11.21 allows for a multiple of applications. In this section we discuss only a small part 
of them but we wish to note that almost all results that have been obtained in the literature 
under restrictions on all jumps of £ can now be strengthened by imposing conditions only on 
positive jumps. This is due to (II. 6h and the fact that on the right-hand side of the identity the 
law of the exponential functionals has been determined by its integral moments which admit 
some simple expressions, see Propositions 14.61 and 14.71 below. 

The factorization allows us to derive some interesting distributional properties. For instance, 
we can show that the random variable Ig is unimodal for a large class of Levy processes. We 
recall that a positive random variable (or its distribution function) is said to be unimodal if there 
exists a S M + , the mode, such that its distribution function F{x) and the function 1 — F(x) 
are convex respectively on (0, a) and (a,+oo). It can be easily shown, see e.g. [28J, that the 
random variable Iy, as defined in Theorem 11.21 is self-decomposable and thus, in particular, 
unimodal. It is natural to ask whether this property is preserved or not for 1^. We emphasize 
that this is not necessarily true even if Ijj- is unimodal itself. Cuculescu and Theodorescu 
|12j provide a criterion for a positive random variable to be multiplicative strongly unimodal 
(for short MSU), that is, its product with any independent unimodal random variable remains 
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unimodal. More precisely, they show that either the random variable has a unique mode at 
and the independent product with any random variable has also an unique mode at or the law 
of the positive random variable is absolutely continuous with a density m having the property 
that the mapping x —¥ logm(e x ) is concave on M. We also point out that it is easily seen that 
the MSU property remains unchanged under rescaling and power transformations and we refer 
to the recent paper [30] for more information about this class of random variables. 

We proceed by recalling that as a general result on the exponential functional Bertoin et al. [31 
Theorem 3.9] have shown that the law of is absolutely continuous with a density which we 
denote throughout by mg. 

In what follows, we show that when £ is a spectrally negative Levy process (i.e. H(dy)I^ y> Qy = 
in (jl.ip and £ is not the negative of a subordinator), we recover the power series representation 
obtained by the second author in [26J for the density of I^. We are now ready to state the first 
consequence of our main factorization. 

Corollary 2.1. Let £ be a spectrally negative Levy process with a negative mean. 

(1) Then, we have the following factorization 

(2.1) I^Ih-xG" 1 , 

where G 7 is a Gamma random variable of parameter 7 > 0, where 7 > satisfies the 
relation ^(7) = 0. Consequently, if Iff- is unimodal then 1^ is unimodal. 

(2) The density function ofl^ has the form 

2.-7-1 f'OO 

(2.2) m c (x) = -—y J e-y/ x y^m H - (y)dy, x > 0, 
where V stands for the Gamma function. In particular, we have 

lim x 7+ mg(x) - 



r(7) 

(3) Moreover, for any 1/x < lim s 



>00 „ ! 



( } d ) ~ r( 7 )r(7 + i) x ^ l) nLi*(fc + 7) x ' 

(4) Finally, for any j3 > 7 + 1, the mapping x <— y x~"m^{x~ 1 ) is completely monotone on 
R + , and, consequently, the law of the random variable LJ 1 is infinitely divisible with a 
decreasing density whenever 7 < 1 . 

Remark 2.2. (1) From pjj Corollary VII. 5] we get that 

lim - [ b ~ I' 1 yU{dy) ifa = ° A y)Ii{dy) < °°' 

s^oo s I +00 otherwise. 

Since we excluded the degenerate cases, we easily check that b — f®^ yll(dy) > 0. 
(2) We point out that in |26j . it is proved that the density extends to a function of a complex 
variable which is analytical on the entire complex plane cut along the negative real axis 
and admits a power series representation for all x > 0. 
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To illustrate the results above, we consider *(s) = — (s — 7 )0_(s),s > 0, with 7 > 0, and 
where for any a € (0, 1), 

(2 .4) _ /<"';-');" 

1 (as + lj 

= / (l-e~ sy ) 1 - 7 r^dy= / (1 - e- sy )w a (y)dy 

Jo K J aT(a + l){ey/ a - l) 2 ~ a 7 

is the Laplace exponent of a subordinator. Observing that the density vr Q (y) of the Levy measure 
of 0_ is decreasing, we readily check that is the Laplace exponent of a spectrally negative 

Levy process. Next, using the identity 1 H - = Gi a , see e.g. [25], we get 
which, after some easy computations, yields, for any x > 0, 

' n=0 

(2-6) = r(7)r (7 + i) ^((a^T + i);-^), 

where 1.F0 stands for the so-called Wright hypergeometric function, see e.g. Section 12.1]. 
Finally, since G\ a is unimodal, we deduce that 1^ is unimodal. Actually, we have a stronger 
result in this case since 1^ is itself MSU being the product of two independent MSU random 
variables, showing in particular that the mapping x h4 ii ? o((a,a7 + l)',e x ) is log-concave on R 
for any a € (0, 1) and 7 > 0. 

We now turn to the second application as an illustration of the situation P+ of Theorem 
11.21 We would like to emphasize that in this case in general we do not require the existence of 
positive exponential moments. We are not aware of general examples that work without such a 
restriction as (14. ip is always crucially used and it is of real help once it is satisfied on a strip. 

Corollary 2.3. Let £ be a Levy process with —00 < E[£i] < and a 2 > 0. Moreover assume 
that 

n(dy)hy>o } = cXe- x ydy, 

where c, A > 0. Then, we have, for any s > — X, 

s 2 

ip+(-s) = 5 + s 2 + k + s + c-——, 

A + s 

where c_ = c/</>_(A) and 5+ > 0. Consequently, the self- decomposable random variable Iy admits 
the following factorization 

(2.7) I r = d + G^ xF^A-^i), 

where < Q\ < A < 62 are the two positive roots of the equation ip+(s) = and B stands 
for a Beta random variable. Then, assuming that 62 — Q\ is not an integer, we have, for any 
1/x < lim^oo \4>-(s)\, 



k + T(X + l)x~ l (J^ E V%- 

m^(x) — 



T(e 1 + i)r(0 2 + 1) \ ^ r(0i + 1) 



x-^Zkiiei + h-x- 1 ) 



where 

oo n 

(2.8) 1^(01 + l;x) = V «„(</>-, 0;)^7 

n=0 

nnH n U ft A - T7 2 r(^-^-n) r(n+0 I+ l) ,■ _ 1 9 

and a n {<p-,V l ) - [[ 3 x _ r(A _ ._ n) nLi 0-^+^) ' ~ ' 

Remark 2.4. TTie assumption a 2 > 0, as ioe// as the restriction on 62 — 61, have been made in 
order to avoid dealing with different cases but they can both be easily removed. The latter will 
affect the series expansion (I5.3p . The computation is easy but lengthy and we leave it out. 

Remark 2.5. The methodology and results we present here can also be extended to the case 
when the Levy measure IL(dy)Ir y>0 \ is a mixture of exponentials as in [9] and [16] but we note 
that here we have no restrictions on the negative jumps whatsoever. 

We now provide an example of Theorem 11.21 in the situation P±. 

Corollary 2.6. For any a £ (0, 1), let us set 

T , x azT(a(— z + 1) + 1) , , „ m 
2.9 = n 1 \ r( -TrcM*), z e iR, 

(1 — z)T(— az + 1) 

where </>+ is as in (|1.2p with /i + 6 V, k + > 0. Then ^> is the Laplace exponent of a Levy process 
£ which drifts to —00. Moreover, the density ofl^ admits the following representation 

T —l/a roo , . 

(2.10) m^x) = — — J g a [(y/x) 1 / a ) mY (y)y 1 / a " 1 dy,x>0, 

where g a is the density of a positive a-stable random variable. Furthermore, i/lim^oo s Q_1 (/>+(— s) 
0, then for all x > 0, 

n=l v 7 

Finally, the positive random variable is MSU if and only if a < 1/2. Hence 1^ is unimodal 
for any a < 1/2. 

Remark 2.7. The fact that l H - is MSU if and only if a < 1/2 is a consequence of the main 
result of [30] . 



Remark 2.8. Note that this is a very special example of the approach of building the Levy 
process from <j>± when fj,± € V . One could construct many examples like this and this allows for 
interesting applications in mathematical finance and insurance, see e.g. 



As a specific instance of the previous result, we may consider the case when 

Tja's + 1) 
T(a'(s + 1) + 1) 

with a' S (0, 1). We easily obtain from the identity (14. 8h below that 

. , _ T(a'm + l-a') 
t[l y j - r(1 _ a/) , m-1,2,..., 

that is Iy = G^" ,. Hence, as the product of independent MSU random variables, 1^ is MSU 
for any a' 6 (0, 1) and a < 1/2. Moreover, using the asymptotic behavior of the ratio of gamma 

7 



functions given in (j5.7|) below, we deduce that for any cJ € (0, 1 — a) we have 

/r,m\ i \ 1 ^T(a'n + 1) 

(2-12) mAx) = — — > — — - (-1) x , 

n=l 

which is valid for any x > 0. 

We end this section by describing another interesting factorization of exponential function- 
als. Indeed, assuming that fi- € V, it is shown in [251 Theorem 1] that there exists a spec- 
trally positive Levy process Y = (Yt)t>o with a negative mean and Laplace exponent given by 
ip,(—s) = —s(j)-(s + 1), s > 0, such that the following factorization of the exponential law 

(2.13) l H - x I" 1 £ G 1 

holds. Hence, combining f)2 . 13|) with (jl.6p . we obtain that 

I C X I^ 1 = Gi X Iy. 

Consequently, we deduce from [29} Theorem 51.6] the following. 

Corollary 2.9. // in one of the settings of Theorem we assume further that £ V , then 
the density of the random variable Ig x I^ 1 , where Iy is taken as defined in ()2. 13[) , is a mixture 
of exponential distributions and in particular it is infinitely divisible and non-increasing on M + . 

Considering as above that = Gf in Corollary 12.11 and 12.31 we deduce from [251 Section 
3.2] that the random variable S~ a x Ig is a mixture of exponential distributions, where S a is a 
positive stable law of index a. 



3. Some results on generalized Ornstein-Uhlenbeck processes 

The results we present here will be central in the development of the proof of our main 
theorem. However, they also have some interesting implications in the study of generalized 
Ornstein-Uhlenbeck processes (for short GOU), and for this reason we state and prove them in 
a separate section. 

We recall that for a given Levy process £ the GOU process U^, is defined, for any t > 0, x > 0, 

by 

rt 

(3.1) Ul{x)=xe^ + e^ e^ s ds. 

Jo 

This family of positive strong Markov processes has been intensively studied by Carmona et 
al. [10] and we refer to [23] for some more recent studies and references. The connection with 
our current problem is explained as follows. From the identity in law &t-s)~)o<s<t = (S,s)s<t, 
we easily deduce that, for any fixed t > 0, 

Uf(x) = xe^ + f e^ 3 ds. 
Jo 

Thus, we have if lim^oo ^ = — oo a.s., that 

uL(x)^k 

and hence the law of le is the unique stationary measure of U^, see [101 Proposition 2.1]. 

In the sequel we use the standard notation G&(IR) (resp. Gfe(M+)) to denote the set of bounded 
and continuous functions on 1R (resp. on M+). Furthermore, we set V' = G^(M), where G^(R) 
is the set of twice continuously differentiable bounded functions which together with its first 
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two derivatives are continuous on 1 = [—00,00]. Then, we recall that, see e.g. [10] for the 
special case when £ is the sum of a Brownian motion and an independent Levy process with 
bounded variation and finite exponential moments and [T7] for the general case, the infinitesimal 
generator L ui of takes the form 

(3.2) L u 'f(x) = L«/ e (lnx) + f'(x), x > 0, 

whenever E[|£i|] < 00 and f e {x) = f(e x ) G Dom(L^), where stands for the infinitesimal 
generator of the Levy process £, considered in the sense of ltd and Neveu (see [201 p. 628-630]). 
Recall in this sense V C Dom(L^) and hence V = {/ : S+ (->• !|/ e G V'} C Dom{L m ). 
In what follows we often appeal to the quantities, defined for x > 0, by 

(3.3) H(x) := / U(dy); U ± (x) := / II ± (dy), 

J\y\>x J y>x 

(3.4) rl(x) := I U(y)dy; Tl ± (x) := [ U ± (y)dy, 

J y>x J y>x 

where U + (dy) = H(dy)l{ y>0 } and U-(dy) = U(—dy)l^ y>0 j. Note that the quantities in (|3.4|) 
are finite when E [|£i|] < 00. Moreover, when E[£i] < 00, (II. ip can be rewritten, for all z G C, 
where it is well defined, as follows 

2 i-oo _ roo _ 

(3.5) tt(z)=E — z 2 + z 2 / e^n+(y)dy + z 2 / e "^n_(y)dy. 

For the proof of our main theorem we need to study the stationary measure of and in 
particular L u ^ in detail. To this end, we introduce the following functional space 

K = {/ : R + M- 1| f e G V; Jirn^ + \f'J(x)\) = 0; / (|£(x)| + |^'(x)|)dx < 00}, 

where / e (x) = /(e x ). 

Proposition 3.1. Let be a GOU process with E[|£i|] < 00. Then K. C Dom{L u<L ). Moreover, 
for any f G K, we have, for all x > 0, 

(3.6) L^/to = £^+E[CiM^) + + / 5'(2/)5 + (ln^)dy+ / «/(y)5_ ( In -) dy, 

where g(x) = xf'(x). Finally, for any function h such that f^iu" 1 A l)\h(y)\dy < 00 and / G /C 
we /iaue 

(3.7) (L u 'f,h) = (g',Ch), 

where (fi,f 2 ) = | °° fi(x)f 2 (x)dx and 
(3.8) 

2 poo / 1 \ /*oo 

£/ l(x ) = ^- x / l (x) + y ( - + EKilJ h(y)dy + J g_ (in |) %)dy + J %■( ln ~) %)^- 

Remark 3.2. There are certain advantages when using the linear operator C instead of the 
generator of the dual GOU. Its integral form allows for minimal conditions on the integrability 
of \h\ and requires no smoothness assumptions on h. Moreover, if h is positive, Laplace and 
Mellin transforms can easily be applied to Ch(x) since the justification of Fubini Theorem is 
straightforward. 
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Proof. Let f E JC then by the very definition of 1C we have that f e € V and from (|3.2|) we get 
that K C Dom(L m ). Next, can be found in [17] but can equivalently be recovered from 
f)3.2f) by simple computations using the expression for L*, which can be found on [H p. 24]. 
To get (|3,7p and (|3.8[) . we recall that g(x) = xf{x) = f' e (lnx) and use (|3.6p combined with a 
formal application of the Fubini Theorem to write 

poo n („.\ poo poo 

(L uS f,h) = / i Mlh{y)dy + - yg , (y)h(y)dy + M[C 1 ] g{y)h{y)dy + 

Jo y 2 Jo Jo 

coo py _ poo poo _ 

g'(v)n_(ki»)dvh(y)dy + / / g'(v)U + (]n-)dvh(y)dy 

JO v JO Jy U 

oo 



pOO L/ A pOO pOO I pOO 

g'(v) -Utldydv + Elh] g'(v) h(y)dydv + — vg'(v)h(v)dv + 

Jv y Jo Jv 2 j 

poo poo _ poo f V = 11 

/ g'(v) U_(hi^)h(y)dydv+ g'{v) / U + (in -)h(y)dydv 

Jo Jv v Jo Jo y 

(3.9) = (g',Ch). 

To justify Fubini Theorem, note that / € /C implies that lim x _>.o g{x) = lim^-j.o fe() nx ) = 0, 
9( x ) = !o d'{v)dv and 

poo r 

/ \g'(v)\dv= / \f'M\dy<C(g)<oo, 
Jo Jr 

(3.10) \g(x)\ + x\g'(x)\ = \f e Qnx)\ + |/^(lnx)| < C(g) < oo, 

where C(g) > 0. Note that f|3. 10|) and the integrability of (1 A y~ 1 )\h(y)\ imply that 



a(y) 



y 



oo ry 

Ky)dy < I / WWdvy-^h^dy <C(g) \ y' l \h{y)\dy < oo, 
o Jo 



and so Fubini Theorem applies to the first term in ()3.9|) . The second term in (j3.9j) remains 
unchanged whereas for the third one we do the same computation noting that only y~ l is not 
present. From (|3.1U|) and the fact that n + (l) + n_(l) < oo since E[|£i|] < oo, we note that for 
the other two terms, we have with the constant C(g) > in (|3,10p . 

j \g'(v)\W_(^ln^jdv = J \xe~ w g'{xe- w )\n_{w)dw 

f' OO pi 

< TT_(1) / \g'(v)\dv + C(g) Tl„(w)dw < oo 

Jo Jo 

[°° \g'(v)\Tt + (]n~)dv = r \xe w g'(xe w )\Ti + {w)dw 

Jx \ XJ Jq 

p oo pi 

< n+(l) / \g'(v)\dv + C(g) H + {w)dw < oo. 

Jo Jo 

Therefore we can apply Fubini Theorem which completes the proof of Proposition 13.11 □ 

The next result is known and can be found in [17] but we include it and sketch its proof for 
sake of completeness and for further discussion. 

Theorem 3.3. Let U$ be a GOU where -oo < E[£i] < 0. Then has a unique stationary 
distribution which is absolutely continuous with density m and satisfies 

(3.11) Cm(x) = for a.e. x > 0. 
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Remark 3.4. Note that due to the discussion in Section^ m = mg, i.e it equals the density 
of the law of Ig . Therefore all the information we gathered for in Section [H is valid here for 
the density of the stationary measure ofU^, i.e. m. 

Remark 3.5. Equation (|3.11j) can be very useful. In this instance it is far easier to be studied 
than an equation coming from the dual process which is standard when stationary distributions 
are discussed. It does not presuppose any smoothness of m but only its existence. Moreover, as 
noted above A3. llf) is amenable to various transforms and difficult issues such as interchanging 
integrals using Fubini Theorem are effortlessly overcome. 

Remark 3.6. It is also interesting to explore other cases when a similar equation to (|3.11[) can 
be obtained. It seems the approach is fairly general but requires special examples to reveal its 
full potential. For example, if L is an infinitesimal generator, M is a differential operator, C is 
an integral operator and it is possible for all f G Co°(K+), i.e. infinitely differ entiable functions 
with compact support, and a stationary density u to write 

(Lf,u) = (Aff,£u) = 

then we can solve the equation in the sense of Schwartz to obtain 

NCu = 0, 

where N is the dual of J\f. If we show that necessarily for probability densities Cu = 0, then we 
can use C to study stationarity. 

Proof. From (|3.7p and the fact that m is the stationary density we get, for all g{x) = xf'(x), 
with / e Cg°(R+) C K, 

(g',£m) = Q. 

Then from Schwartz theory of distributions we get Cm{x) = Clnx+D a.e.. Integrating (13.8P and 
the right-hand side of the latter from 1 to z, multiply the resulting identity by z" 1 , subsequently 
letting z —¥ oo and using the fact that m is a probability density we can show that necessarily 
C = D = 0. The latter requires some efforts but they are mainly technical. □ 

Theorem 3.7. Letm be a probability density function such that f^°m(y)y~ 1 dy < oo and ([3. lip 
holds for m then 

(3.12) m{x) = m{x) a.e., 

where m is the density of the stationary measure of . 

Remark 3.8. This result is very important in our studies. The fact that we have uniqueness 
on a large class of probability measures allows us by checking that (|3.1ip holds to pin down the 
density of the stationary measure ofU^ which is of course the density ofl^. The requirement that 
Jq 00 m(y)y~ 1 dy < oo is in fact no restriction whatsoever since the existence of a first negative 
moment of Lj is known from the literature, see [3] . 

Remark 3.9. Also it is well known that if L u is the generator of the dual Markov process then 
L u m = does not necessarily have a unique solution when L u is a non-local operator. Moreover 
one needs assumptions on the smoothness of m so as to apply L u . Using C circumvents this 
problem. 

Proof. Let (Pt)t>o be the semigroup of the GOU U*, that is, for any / € Cb(M+), 



Ptfix) = E [/ (Ul{x) 
li 



x > 0, t > 0. 



If (|3.1ip holds for some probability density m then ()3.7p is valid, i.e. for all / G fC, 

(L u 'f,m) = (g',Cm)=0. 

Assume for a moment that 

(3.13) P S K C K, for all s > 0, 

and, there exists a constant C(/, £) > such that, for all s <t, 



(3.14) 



L uS P s f(x) <C(f,t)(x- 1 Al). 



Then integrating out with respect to m(x) the standard equation 

P t f(x) = f(x)+ f L u 'PJ(x)ds, 
Jo 



we get, for all / G /C, 



Ptf(x)m(x)dx 



f{x)m(x)dx. 



Since Co°(R + ) C /C and Cq°(R + ) is separating for Co(R+), the last identity shows that m is 
a density of a stationary measure. Thus by uniqueness of the stationary measure we conclude 
([3TT2]) . Let us prove (^T3l) and (I3T4"]) . For / G /C write 



['(< 




= E 

















g s (x) :=P s f(x)=E 

Put 5 S (x) = g s (e x ) = (g s ) e (x) . Note that since / G K. and < e x+ ^ s < e x+ ^ s + J Q S e^ v dv we have 
the following bound 



(3.15) 



e x+ ^f (e x+ t° + J'^dvj 


+ 











<C(f) 



which holds uniformly in x G R and s > 0. In view of (|3.15p the dominated convergence theorem 
gives 



g' s (x)=E 



^+6 y/ / g x+£ 



+ 



m^{\~g' s {x)\,\~£{x)\} <C{f). 



+ E 



dv 



(3.16) 

Clearly then from ()3.15p . ()3.16p . the dominated convergence theorem and the fact that / G K 
which implies the existence of lim a ._ KX) f"(x) = b, we have 



lim g£(x) = E 



lim 

x— >oo 



Similarly, we s how t h at lim ^oo g' s (x) = lirn^oo f' e (x) and trivially lim a .- + ± o g s (x) = lim^ioo f e (x). 
Finally using (|3.15p . ([3.16p . / G /C, the dominated convergence theorem and the fact that for all 
s > almost surely J* Q S e^ v dv > 0, we conclude that 



lim \g' s (x)\+\g'^x)\ < 2E 



lim 

x— >— oo 



D x+is f> I e x+i 



+ 



^ v dv 



+ 



2(a-+6) J" I p x+t 



+ 



which together with the limits above confirms that g s G V' and proves that 



Urn _\g' s (x)\ + \g'^x)\ =0. 
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Finally since / £ K and (j3. 15|) . we check that 



and 



\g'M\dy<M 



\~g'l{y)\dy < E 



Io eivdv 



\f(u)\du 



< 



\f'(u)\du = / \f' e {u)\du<CU), 



u- e^dv) \ f"(u)\du 



< 



u\f"(u)\du 



< 2 f \f e (\nx)\ + |#(]n*)| — = 2 f (\f e (y)\ + \f^(y)\)dy < C(f), 

where C(f) is chosen to be the largest constant in all the inequalities above and we have used 
the trivial inequality u 2 \f"(u)\ < |/^(lnu)| + |/^(lntt)j. Thus using all the information above 
we conclude that g s = P s f & K- and (|3,13p holds. Next we consider (|3.14p keeping in mind that 
all estimates on g s we used to show that g s € /C are uniform in s and x. We use (|3.6|) with 
g(x) = xg' s (x) = g' s (lnx), the bounds on g s and its derivatives to get 

^ + E[Zi]g(x) + yx</(x)| < C(f)x^ + C(/)|E[fr]| + C(/)y < C(f, a, E[fr])(l A x^ 1 ). 
Moreover, as in the proof of Proposition 13.11 we can estimate 



g {v)Ii- (In -)dv 



+ 



g'(v)n + (\n-)dv 

x 



< 



(n_(i) + n+(i)) J™ \ g '(s)\ds + c(f) Qf 1 U-(s)ds + J n H 



(s)ds 



n_(i) + n+(i) 



and therefore (|3.14p holds since 



g"(y)\dy + C(f)[ / n_(a)ds + / n + (s)ds) <C 



L u 'g s (x) = ^+n^]g(x) + ^-xg'(x)+ \ [" g'(v)IL_ (in ~ )dv + / ff»n + (ln 



This concludes the proof. □ 
Theorem 3.10. Let {^ n ^) n >i be a sequence of Levy processes with negative means such that 

lim e (n) = e, 

where ^ is a Levy process with E[£i] < 0. Moreover, if for each n > 1, stands for the law 

of the stationary measure of the GOU process U defined, for any t > 0, x > 0, fry 



5 («) ,(») 



and the sequence (m^) n >i is tight then {m^) n >i converges weakly to m^, which is the unique 
stationary measure of the process , i.e. 



(3.17) 



lim mS n ^ 



m 



(o) 



13 



Proof. Without loss of generality we assume using Skorohod-Dudley theorem, see Theorem 3.30 
in Chapter 3 in [15], that the convergence e (n) -> £ holds a.s. in the Skorohod space £>((0, oo)). 
Due to the stationarity properties of mf- n \ for each t > 0, we have, for any / £ C&(R+), 

(/,mW) = (P t (n) /,m (n) ) = (p^f-PtfM^) + (PtfM n) ) , 

where P+ and Pf are the semigroups of and t/ £ . For any x > 0, 



(p/ re) /-p/, 



(n) 



< 211/IUmW (a-, oo) + sup |P, (n) / (y) - P t / (y) 

f/<rr 



(3.18) 



< 2||/|| m (n) (z, oo) + E 



sup 

y<x 



f (y)) - / (V| (y)) 



Taking into account that (m( n )) n >i is tight we may fix 5 > and find x > big enough such 
that 

supm'"'(j;, oo) < 5. 

Also since / € Cfc(K+) then / is uniformly continuous on K + . Therefore, to show that 



lim E 

n— >oo 



sup 

y<x 



f{ut\y))-f(uHy)) 



0. 



due to the dominated convergence theorem all we need to show is that 



(3.19) 



lim sup \uf n \y) - Uf(y)\ = 0. 
n ^°°y<x 



From the definition of U^ (n) and IJt, we obtain that, for y < 



e(n) 



US (y)-UHy) 



< x 


An) , 


+ 


An) * 















Jo j o 



e ?s - e K "ds 



Since £( n ) ^4' £ in the Skorohod topology and 

P ({3n > 1 : C t (n) - e£ } > 0} n {6 - 6- > 0}) = 

the first term on the right-hand side of the last expression converges a.s. to zero asn-> oo. The 
a.s. convergence in the Skorohod space implies the existence of changes of times (A n ) n >i such 
that, for each n > 1, A n (0) = 0, X n (t) = t, the mapping s t— >■ X n (s) is increasing and continuous 
on [0, £], and 



(3.20) 



lim sup |A n (s) — s\ = lim sup |A„ (s) — s| = 

n->oo s<t n-)-oo s<t 



(3.21) 
Hence, 



lim sup 



*•(») _ £ 



lim sup 



0. 



e ^ — e ^ s as 



< 



An) 



li 



+ 



g ^An 1 ^) _ g 



The first term on the right-hand side clearly goes to zero due to (|3.2ip whereas (|3.20p implies 
that the second term goes to zero a.s. due to the dominated convergence theorem and the fact 
that pathwise, for s < t, 

' > 



lim sup 

n— >-oo 

only on the set of jumps of £ and this set has a zero Lebesgue measure. Thus we conclude that 

rt 

e — e 



lim 

n— >oo 







0. 



Similarly we observe that 



(3.22) 



lim 

n— >oa 



(n) 



e ?s ds < lim i 



;(n) 



where the last identity follows from 



sup 



(-ei n) ) 



< sup 



+ sup 

s<t 



AnOO 



SUp \Cs\+ SUp 



A„ (s) 



and an application of (|3.2ip . Therefore, f|3. 19|) holds and 



lim sup 

n ^°° y<x 



f(ui n) (y))-f(uHy)) 



o. 



The dominated convergence theorem then easily gives that the right-hand side of (I3.18P goes to 
zero and hence 



lim sup 



P t {n) f -P t f,m^) < 2| |/| loo sup (x, oo) < 2 

n>l 



As (5 > is arbitrary we show that 



lim 

n— >oo 



P{ n) f-P t f,m^ 



0. 



Since (m.( n )) n >i is tight we choose a subsequence (m^ nfc ^)fc>i such that lim^oo m^ nk ^ = v with 
f a probability measure. Then, for each t > 0, 

C/»= lim (/m^) = lim ( P t {nk) f, m {nfc) 



lim (P t /,m (nfc) 



Therefore f is a stationary measure for U^. But since m^ ' is the unique stationary measure we 
conclude that 

lim mS n ^=v = . 



This translates to the proof of (13.171) . 



□ 



4. Proof of Theorem Q 

We start the proof by collecting some useful properties in two trivial lemmas. The first one 
discusses the properties of VP. 

Lemma 4.1. \29\ Theorem 25.17] The function defined in (|3.5p . is always well-defined on 
iR. Moreover, 4' is analytic on the strip {z G C; — a_ < $t(z) < a+}, where a_,a + > if and 
only ifE [e(- - +e )fr] < oo and E [e( a +- e )? 1 ] < oo for all < e < a_ A a+. 

The second lemma concerns the properties of (f>± and is easily obtained using Lemma 14.11 
(|1.4p together with the analytical extension and the fact that subordinators have all negative 
exponential moments. 
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Lemma 4.2. Let £ be a Levy process with E[£l] < oo. T/ien 0+ is always analytic on the strip 
{z € C; 9ft(z) < 0} and is well-defined on iR. Moreover <j) + is analytic on {z G C; < a+}, /or 
a + > 0, if and only i/E [V a+_e ^ 1 ] < oo, for some e > 0. Similarly 0_ is always analytic on the 
strip {z S C; > 0} and is well-defined on iR and is analytic on {z € C; < — a_}, 
/or a_ > 0, if and only if E[e^ _a_+<E ^ 1 ] < oo, for some e > 0. Finally, the Wiener-Hopf 
factorization (jl.4p /io/ds on i/ie intersection of the strips where ^+ and 0_ are well-defined. 

4.1. Proof in the case E+. We recall that in this part we assume, in particular, that £ is a 
Levy process with a finite negative mean and that there exists a + > such that |^(^)| < oo for 
any < 3?(z) < a + . Next, we write 9* = max(0,a+), where 9 = inf{s > 0; = 0} (with the 
convention that inf = +oo). We also recall from |10j . see also [21j . that the Mellin transform 
of i£ defined by 

/•oo 

■M-mAz) = I x z ~ l m^{x)dx 
Jo 

satisfies, for any < < 9*, the following functional equation 

(4.1) M mi (z + l) = -^M ms (z). 
We proceed by proving the following easy result. 

Lemma 4.3. If /i + € V then there exists a spectrally positive Levy process Y with Laplace 
exponent ip + (—s) = —s<j>+(—s), s > 0, and a negative finite mean —<j) + (0). Moreover, if £ has a 

negative finite mean then E (I#- x Iy)" 1 = —(p+(0)(p'_(0 + ) < +oo. 

Proof. The first claim follows readily from jTJ Theorem VII.4(ii)] and by observing that ^+(0 - ) = 
0+(O). From (|4.8j) we get that E [iy 1 ] = k + . Next, since — oo < E[£i] < 0, using the dual version 
of [131 Corollary 4.4.4(iv)], we get that -oo < E[i/f ] < and thus 0'_(O + ) < oo. From the 
functional equation (|4.ip . we easily deduce that E [i^ 1 -] = — 4>'_(0 + ) which completes the proof 
since the two random variables are independent. □ 

Lemma 4.4. Assume that £ has a finite negative mean and condition E + holds. Let rj be a 
positive random variable with density k(x), such that E^" 1 ] < oo and E[n 5 ] < oo, for some 
9* > 5 > 0. Then, for any z such that ^ft(z) S (0, 5), 

f°° , v&fV) 1 

(4.2) M Ck (z)= / x z ~ 1 CK(x)dx = A I L M K (z + l) + -M K {z) 

Jo z 2 * 

and if Mck(z) = 0, for < a < $l(z) < b < 5, then Ck(x) = a.e.. 

Furthermore the law of the positive random variable Iy x I H - , as defined in Theorem \l.2\ is 
absolutely continuous with a density, denoted by m, which satisfies 

(4.3) Cm(x) = for a.e. x > 0. 

Remark 4.5. Note that the proof of this lemma shows that we have uniqueness for the probability 
measures with first negative moment that satisfy (|4,ip . This is a rather indirect approach and 
seems to be more general than the verification approach of [16] . see Proposition 2, where precise 
knowledge on the rate of decay of the Mellin transform M m( (z) is needed. In general such an 
estimate on the decay seems impossible to obtain. 

Proof. We start by proving (|4.2p . Note that since f °° y~ l n{y)dy < oo, we can use Proposition 

ED to get 

(4-4) 

C K {x) = ^-xk{x)+ l°° ^Mdy+E^] I" n{y)dy+ /"°°IL (in*) K(y)dy+ f 5+ fin-) n{y)dy. 

1 J x y Jx Jx v x/ Jo v yj 
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As k is a density, one can use Fubini Theorem to get, after some easy computations, that for 
any e < < 5 < 9*, with < e < 5, 



oo 

Mck(z) = I x z ~ 1 Ck(x)<1x 
Jo 

= M K {z + 1) + + j^M^e^dy + r! + (y)e 2 ^y) + 



^fzl 1 



Let -M£ K (z) = for e < < 5. We show using that all terms in (|4.4|) are positive except the 
negative one due to E[£i] < that, with u = $l(z), 

rV" 1 \£ K (x)\dx < M K {u + 1) -2^-) +-M K {u) <oo. 

Jo \ u z u J u 

Given the absolute integrability of x z ~ l Ck(x) along imaginary lines determined by e < ^t(z) < 5 
we can apply the Mellin inversion theorem to the identity Mck{z) = to get £n(x) = a.e., 
see Theorem 6 in Section 6 in [8]. 

Next it is plain that the law of Iy x I H ~ is absolutely continuous since, for any x > 0, 

f°° / x\ 

(4.5) m(x) = J m Y y-jy~ 1 m H ~(y)dy. 

Furthermore from the Wiener-Hopf factorization (jl,4p and the definition of we have that 

— z —z —z 



which is valid, for any < < 9*. Thus, we deduce from the functional equation (|4.ip and 
the independency of Y and H~ that, for any < 5R(z) < 9*, 

(4.6) M m (z + l) = --^-M m (z). 

Next, since from Lemma [43} we have that y~ 1 m(y)dy < oo, we can use Proposition 13.11 and 
thus (|4.4p and subsequently (|4.2p are valid for m. Moreover due to the representation (|3.5p of 
* and relation (|4.6|) we have that for any e < $l(z) < 9* with < e < 9*/ 4, 

Mcm(z) = ^AU> + 1) + -M m (z) = 

z z z 

and we conclude that £fn{x) = a.e. 

□ 

We are now ready to complete the proof of Theorem 11.21 in the case E+. Indeed, since m^, 
the density of It, is the density of the stationary measure of we have that is also solution 
to (|3.1ip . Combining Lemma 14.41 with the uniqueness argument of Theorem 13.71 we conclude 
that the factorization (11.61) holds. 



4.2. Proof of the two other cases : P+ and P±. We start by providing some results which 
will be used several times throughout this part. 
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Proposition 4.6 (Carmona et al. [_0]). Let H be the negative of a (possibly killed) subordinator 
with Laplace exponent <f>, then the law of Ih is determined by its positive entire moments as 
follows 

(4-7) Epg] = ^ + * } m = l,2,■■■ 

llfc=l (-n k )) 

Proposition 4.7 (Bertoin and Yor [5]). Let Y be an unkilled spectrally positive Levy process 
with a negative mean and Laplace exponent ip + , then the law o/l/Iy is determined by its positive 
entire moments as follows 

(4.8) E[ly m ] = E[-y 1 ] n ^ V ' +( ~ A;) ,m = l,2,..., 

T(m) 

with the convention that the right-hand side is E[— Y{\ when m = 1. 

In order to get (|1 .6j) in the case when £ does not have some finite positive exponential moments, 
we will develop some approximation techniques. However, the exponential functional is not 
continuous in the Skorohod topology and therefore we have to find some criteria in order to 
secure the weak convergence of sequences of exponential functionals. This is the aim of the next 
result. 

Lemma 4.8. Let (^ n ^) n >i be a sequence of Levy processes with negative means such that 

lim C (n) = . 

n— ¥00 

where £ is a Levy process with E[£i] < 0. Let us assume further that at least one of the following 
conditions holds: 

(a) for each n > 1, andt; are unkilled spectrally positive Levy processes such that limn._s.oo E[£j n ^] 
E[.i], 

(b) for each n > 1, £W and £ are the negative of unkilled subordinator s, 

(c) the sequence (m^ n )) n >i is tight, where m^ n ) is the law o/I^(„). 

Then, in all cases, we have 

(4.9) lim L (n) = L. 

n— >oo 5 

Proof. To prove fj4.9j) in the case (a), we simply observe that writing ip^ for the Laplace 
exponent of ^\ we have, by Levy continuity Theorem, see e.g. [14} Theorem XIII. 1.2], that for 
all s > 0, ip+ (— s) — > tp + (—s) as n — > oo. Next, putting Mm for the sequence of negative entire 
moments of I^(n), we easily deduce, from (|4.8p for all m = 1,2... , that lim ra _ >00 Mm — M m 
where M m is the sequence of negative entire moments of i£. These random variables being 
moment determinate, see Proposition 14.71 we conclude (a) by invoking [141 Examples (b) p. 269]. 
The second case follows by applying a similar line of reasoning to the expression (14. 7ft . Finally, 
the case (c) is a straightforward consequence of (13.17|) of Theorem 13.101 □ 

Before stating our next result, we need to introduce the following notation. Let us first recall 
that the reflected processes {K£ = sup 0<s<t £ s — £t) t>0 and [Rt = — info< s <t Cs) t>0 are Feller 
processes in [0, oo) which possess local times L^ = (L^) t > at the level 0. The ascending and 
descending ladder times, Z ± = (Z ± (t))t>o, are defined as the right-continuous inverses of L ± , i.e. 
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for any t > 0, l^(t) = inf{s > 0; Lf > t} and the ladder height processes H + = (H + (t)) t > and 
-H- = (-H-(t)) t > by 

H + (t) = Ci+(t) = SU P £s , whenever l + (t) < oo , 

0<s</+(t) 

—H~(t) = Ci-(t) = inf £s , whenever l~(t) < oo . 

; 0<s<l-(t) 

Here, we use the convention inf{0} = oo and H + (t) = oo when L+ < t and —H~(t) = — oo 

when < t. From [I3j p. 27], we have, for a, j3 > 0, 

(4.10) 

-aZ+(l)-/3#+(l) 



logE 



-k(a,P) = -k + - V+ a-5 + P-j^ J [l-e-^ + ^) f i + (dy l ,dy 2 ), 

where n+ is the drift of the subordinator l + and /i+(dyi, dy 2 ) is the Levy measure of the bivariate 

subordinator (l + ,H + ). Similarly, for a, ft > 0, 

(4.H) 

logE J e -(«»-(i)-^-(i))J = -k*(a,p) = -r]_a-6-P- J™ J™ (l - e~ (ayi+ ^ ^{d Vl , dy 2 ), 

where n_ is the drift of the subordinator l~ and /U_(dyi, dy 2 ) is the Levy measure of the bivariate 
subordinator (l~ , —H~). 

Lemma 4.9. Let £ be a Levy process with triplet (a, a, LT) and Laplace exponent tp. Let, for any 
n > 1, £( n ) be the Levy process with Laplace exponent denoted by ip^ and triplet (a, a, H^) 
such that Il( n ) = IT on R_ and on R + 

U^(dy) = h^(y)U(dy), 

where for all y > 0, < h^ n \y) "f 1 as n — >■ oo and uniformly for n > 1 we have that for some 
C > 0, limsup y _ > . y _1 (l - /i n (y)) < C. T/ien, 

(4.12) lim e (n) = f, 
and for all a > 0, f3 > 0, we have, as n — > oo, 

(4.13) fc( n ')(a,/5) -> fc(a,/3), 

A;i n) (a,/3) -»• fe»(a,/3), 

where fc( n )(a,/3) and ki (a,(3) stand for the bivariate Laplace exponents of the ladder processes 
of normalized such that k( n \l,0) = k* (1,0) = 1. k(a,(3) and k*(a,/3) stand for the 
bivariate Laplace exponents of the ladder processes of£, normalized such that k(l,0) = fc*(l,0) = 
1. 

Remark 4.10. Denote by (^resp. [l^y—H^yj the bivariate ascending (resp. de- 

scending) ladder processes of^ n ^ and (l + , H + ) {resp. —H~)) the bivariate ascending (resp. de- 
scending) ladder processes of£, then from the Levy continuity Theorem we deduce that asn — > oo, 

0t„).<)4(* + ,# + ), 

where in the convergence the sequence of killing rates of the ladder height processes also converge 
to the killing rate of the limiting process. 
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Proof. For the sake of completeness and also to include both the compound Poisson case and 
the cases when a = and/or (3 = 0, we must improve the proof of Lemma 3.4.2 in [31]. Next, 
since H^ n \dx) A II (dec), where A stands for the vague convergence, we get (|4.12j) from e.g. [I~5l 
Theorem 13.14(i)]. We note the identity 

(4.i5) e = f (n) +l (n) , 

where £^ n ) is a subordinator with Levy measure n(«)(dy) = (1 - h n (y))U(dy) and no drift, since 
1 — h n (y) = 0(y) at zero. Then, when £ is a compound Poisson process we have that £( n ) is a 
compound Poisson process and, for all t > 0, 

£W = o) = P U t = 0, t < f( n) ) + P = 0, t > T (n) 



where f(") = inf{s > 0; || n) > 0}. Since for all y > 0, /i (ri) (y) t 1, then P(t > fW) ^ as 
rt — )■ oo and 

P G dy) I { ,> 0} 4P(& € dy)I {y > 0} . 

When £ is not a compound Poisson process, the law of £W does not charge {0} and thus as 

n — > oo 



^ V ' t; G dyj I { , >0} 4P(&€ dy)I { , >0} . 
Henceforth, from the expression 

y dt J (e"* - e - Q '-^) t^P^ G dy)J 

which holds for any a > and /3 > 0, see e.g. [TJ Corollary VI. 2. 10], we deduce easily that for 
both cases 

(4.17) lim k^ n \a,j3) = k(a,P). 
Moreover, we can write 

(4.18) k^(a,/3) = k( n \0,0) + ~k^(a,f3), 

where k^ are the Laplace exponents of unkilled bivariate subordinators, see [131 p. 27]. Note 
from ([USD that 

fcW(0, 0) = exp (- (1 - e"<) P > 0) *£\ . 

Next from (jiTToT) and the fact that ^ is a subordinator, we have that P > o) < P (& > 0) 

and appealing to the monotone convergence theorem we get that fc( n )(0,0) 1 jfc(0,0). Hence 
we deduce from (|4.17j) and (|4.18j) that for any a, P > 0, U n \a,/3) -»■ fc(a,/8) where fc(a,/3) = 
k(a, (3) — /c(0,0). From the Levy continuity theorem, we have, writing J for the un- 



killed versions of the ascending bivariate ladder processes, that ( Z^, — > ( / + , if + ) , where 



l + , H + j stands also for the unkilled version of y, + ,H + j . These probability distributions being 

proper, we have that for all a, ft G M, k^ n \ia, 0) — > k(ia,i/3), see [HI Theorem XV. 3. 2]. Hence 
k( n \0,i(3) — > k(0,i/3) for all /3 G M which completes the proof for the ascending ladder height 
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processes. The proof of the convergence of the Laplace exponent of the bivariate descending 
ladder process follows readily from the identities 

(n) (i/3)-a = -k( n \a,-iP)ki n) {a,iP) 
ip(i/3) — a = —k(a,—i/3)k*(a,i/3) 

and the convergence of ip^ to ip and k^ to k. □ 

4.2.1. The case P+. We first consider the case when £ satisfies both the conditions P+ and 
> — oo. We start by showing that the condition P+ implies that /i+ £ V . To this end, we 
shall need the so-called equation amicale inversee derived by Vigon, for all x > 0, 

/•oo 

(4.19) fi+(x)= / U + (x + y)U-(dy), 

Jo 

where W_ is the renewal measure corresponding to the subordinator H~ , see e.g. [131 Theorem 
5.16]. 

Lemma 4.11. Let us assume that H + (x) has a non-positive derivative vr + (x) defined for all 
x > and such that —ir + (x) is non-increasing. Then £l+(x) is differentiate with derivative u(x) 
such that —u(x) is non-increasing. 

Proof. Fix x > and choose < h < x/3. Then we have the trivial bound using the non- 
increasing property of —ir + (x) and the description (|4.19p of fi+(x) 

\fi+(x±h)-ji+(x)\ < r \Tl + (x + y±h)-Ti + (x + y)\ 



h Jo h 

POO 

< / (-ir+(x + y-h))U-(dy) 



We show now that the last expression is finite. Note that 

(y + v)) U-(d y ) ^ E (y + n ) ^ n + x ) - • 



From the trivial inequality U^{n + 1) — U-{ri) < £/_ (1), see [18, Chapter 2, p. 11], and since 
— 7r + (x) is the non-increasing density of 11+ (a;), we have with C = U—(l) > 0, 



f°° (Ix \ 
J -n+i-j + yJU-idy) < 






-7T+ | 




(2x 








(2x 







n>l 



3 7 V 3 



Therefore, for all x > 0, the dominated convergence applies and gives 



/■oo 

u(x) = / ir + {x + y)U-{dy). 
Jo 



10 

As — 7T + (x) is non-increasing we deduce that —u(x) is non-increasing as well. □ 
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In the case P+, in comparison to the case E + , we have that £ does not necessarily have some 
positive exponential moments. To circumvent this difficulty we introduce the sequence of Levy 
processes £^ n ) obtained from £ by the following construction: we keep the negative jumps intact 
and we discard some of the positive ones. More precisely, we thin the positive jumps of £ to get 
a Levy process with whose density has the form 

(4.20) „%\x) = ir + (x) (l {0<x < 1} + e-"- 1 ^- 1 )!^) . 

Clearly, —ir+\x) is non-increasing and E e s ^i ' < oo, for s G (0,n _1 ), see (|4.20p . Moreover, 
since we have only thinned the positive jumps and pointwise lim n _ 5 . 00 71+ (x) = tt + (x) , see (|4.20p , 

(4.21) lim C H =' £ 

n—too 

almost surely in the Skorohod space 2?(0,oo). Finally, since -oo < E |fj n) l < E [ft] < and 

— ir^\x) is non-increasing then Lemma 14.111 applies and we deduce that the Levy measure of 
the ascending ladder height process of f M has a negative density whose absolute value is non- 
increasing in x. Then since, for each n > 1, £( n ) has some finite positive exponential moments, 
we have that 

(4.22) L( n ) = 1 H - X Iy(n). 

(n) 

(n) 

Since we thinned the positive jumps of ft for all t > 0, Q < ft and the monotone convergence 
theorem together with (|4,2ip imply that 

(4.23) lim L ( „) =' h. 

n— >oo 5 

By the choice of the approximating sequence we can first use Lemma 14.91 to get 

(4.24) lim H7 £ H~ 

n— too 

and then Lemma 13.101 (b) to obtain that 

(4.25) lim 1^ =I H -. 

Again from Lemma [4.9l we deduce that k^(0, —s) — > k(0, —s), for all s > 0, and lim n _ i . 00 E[Yj^] = 
— linin.j.oo k^ (0, 0) = E[Yi], so we can apply Lemma f3. 101 (a) to get that 

lim Iy(„) = Iy, 

n— >oo 

which completes the proof in this case. □ 



4.2.2. The case P±. First from the philanthropy theory developed by Vigon [31], we know that 
the conditions fi + € V and /i_ 6 V ensure the existence of a Levy process £ with ladder 
processes H + and H~ and such that the Wiener-Hopf factorization (|1.4p holds on zR. Since we 
also assume that k + > 0, this Levy process necessarily drifts to — oo. Next let us introduce the 
Laplace exponents 

(4.26) $\z) = S + z+ [ (e zx -l)^\dx)-k + , 

J(0,oo) 

(4.27) ^\z) = -5-z- [ (1 -e~^)^L n) (dx), 

J(0,oo) 
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where we set fjf\dx) = e x / p fi + (dx), p > 0, and fj,_(dx) = e x l n p, + (dx), n > 0. Plainly, for 
any p > 0, n > 0, € "P and ^i' G T 3 , hence there exists a Levy process £( p ' n ) with Laplace 
exponent \&( p ' n ) satisfying 

(4.28) ^( p ' n )(z) = -<t^\z)(lt\s), 

which is easily seen to be analytic on the strip —1/n < < 1/p. Moreover, from [13\ 

Corollary 4.4.4], we have E[^ p,n) ] = -k+ (/* xe x l n [i+(dx) +5—), which is clearly finite and 
negative. Hence the conditions E + are satisfied and we have, with the obvious notation, that 

If(p,n) = ijj- X Iy(p) 

(n) 

where for any p > 0, is a spectrally positive Levy process with Laplace exponent ipf (— s) = 

—scj)+\—s), s > 0. Let us first deal with the case n — > oo. Since 0_ (s) — > </>_(s), for all s > 0, 
we have that 



lim = £T 



and from Lemma 13.101 (b) we get that 



lim 1^- = l H - . 

n—>oo (n) 

Thus, we deduce that, for any fixed p > 0, the sequence {l^( P , n )) n >i is tight. Moreover, for any 

fixed p > 0, we also have 4 £ (p) , as n — > oo, where £^ has a Laplace exponent given 
by 

(4.29) ^( p )(z) = -<$(z)<f>-{z). 

Indeed this is true by the philanthropy theory. Then from Lemma 13.101 (c), we have that 

lim L( P ,n) = L( P ) = Ih- x Iy(p) , 

which provides a proof of the statement in the case P± together with the existence of some finite 
positive exponential moments. Next, as p — > oo, 4>+\s) — > <A+(s), for all s > 0, and we have 
that 

lim = Y, 

p— >oo 

where Y is a spectrally positive Levy process with Laplace exponent ip+(— s) = — s0+(— s). As 
EfY^] = 4>±\o) = we can use Lemma [3.101 (a) to get 

lim Iy( p ) = Iy. 

p— S-oo 

As above, we conclude from Lemma f3. 101 (c) that 

lim L (p) = I € = 1^- x I y , 

p— >oo s 

which completes the proof of the theorem. □ 
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5. Proof of the corollaries 

5.1. Corollary 12.11 First, since £ is spectrally negative and has a negative mean, it is well 
known that the function ^ admits an analytical extension on the right-half plane which is 
convex on 1R + drifting to oo, with ^'(0 + ) < 0, and thus there exists 7 > such that ^(7) = 0. 
Moreover, the Wiener-Hopf factorization for spectrally negative Levy processes boils down to 

s — 7 

It is not difficult to check that with 4>+(s) = s — 7 and <j>-(s) = —73^, we have //_, € V. 
Observing that ip+(s) = s 2 — 7s is the Laplace exponent of a scaled Brownian motion with a 
negative drift 7, it is well-known, see e.g. [32], that 

l Y = G-\ 

The factorization follows then from Theorem Un considered under the condition P-|-. Since the 
random variable G~ l is MSU, see [T2J, we have that if I H - is unimodal then 1^ is unimodal, 
which completes the proof of (1). Next, (2) follows easily from the identity 

1 f°° 

(5.1) m^x) = ^-yx" 7 " 1 ^ e- y/x y^m H -(y)dy 

combined with an argument of monotone convergence. 

Further, we recall that Chazal et al. Theorem 4.1] showed, that for any (3 > 0, (j)p(s) = 
-^fp4>-( s + P) is a ls° the Laplace exponent of a negative of a subordinator and with the obvious 
notation 

(5.2) rn H -(x) = ^ ± - L , x > 0. 

E[I^_] 

Then, assuming that 1/x < hin^^oo ^(u)/u, we have, from (|4.7p . (|5.ip and (|5.2j) . 

-I OO ^ /*00 

\1) „—n ^" Jo 



n=0 

3 — n 1 



f [JT 1 00 

M a -7-iyf-D 

r w ^0 n*=i-^-(*+7) 

E ^-] fv ir r(n + 7 + i) 

r( 7 )r( 7 + i) ^ ; IK=i-^-(* + 7) 

x -7-i fv ir r(n + 7 + l) 

r(7)r(7 + i) ^ j IK=i*(* + 7) ' 



where we used an argument of dominated convergence and the identity —k<f)_(k + j) = ^/(k + r y). 
Next, again from (|5.ip . we deduce that 



x ^mAx l ) = — ^— x 1+1 13 I e xy y 1 m H -(y)dy 
r (7J 7o 

from where we easily see that, for any /3 > 7 + 1, the mapping x 1— >• x _ ^m^(x _1 ) is completely 
monotone as the product of two Laplace transforms of positive measures. The proof of the 
Corollary is completed by invoking [29^ Theorem 51.6] and noting that LJ 1 has a density given 
by x~ 2 m^(x ), i.e. with (3 = 2. 

24 



5.2. Corollary 12.31 We first observe from the equation (|4.19j) that, in this case, 

/•oo 

Jx+{x) = ce~ Xx / e- Xy U-{dy) 



C-e~ Xx , 



where the last identity follows from [13] and we have set c_ = ^ ^ . From (|1.5p . we deduce 
that Y is a spectrally positive Levy process with Laplace exponent given, for any s < A, by 



s 2 

tp+(s) = 5+s 2 — k + s + c 



A-s 



A-s 

(-5+s 2 - (<5+A + £;+ + c_)s - k + \) , 



where 6+ > since a > 0, see |13|. Corollary 4.4.4]. Thus, using the continuity and convexity 
of ip + on (— oo, A) and on (A, oo), studying its asymptotic behavior on these intervals and the 
identity ip'+(0) = — k + < 0, we easily show that the equation ■0+( s ) = nas 3 roots which are 
real, one is obviously and the two others 8\ and 02 are such that < 9\ < A < 02- Thus, 

M-s) = -T-— — (s + 8i) (s + 6 2 ) , s > —X. 
A + s 

Next, from (|4.8p . we have, with C = k + r^+i^r^+i) an d for to = 2, ... , that 

jgrj-mi = C5 m-I r ( m + gl)£(m + fl 2 ) 

1 y J + T(m + A) 

from where we easily deduce (|2.7p by moments identification. Note that a simple computation 
gives that 0\02 = 5+Xk+ securing that the distribution of Iy is proper. Next, the random 
variable Iy 1 being moment determinate, we have, for ^ft(z) < 6\ + 1, that 

TZ _n ™_ z r(-* + e 1 + i)r(-z + e 2 + i) 



r(-z + A + i) 

Applying the inverse Mellin transform, see e.g. |22[ Section 3.4.2], we get 
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(5.3) my (^) = cE x " e<_1 ^(- x_1 )' :c>0 ' 



i=l 

where 2j(x) = b n,i^, & M = ^x-e^-n) and V 2 = ^x-e^-n) ■ The P roof of the Corollary 

is completed by following a line of reasoning similar to the proof of Corollary 12.11 

5.3. Corollary 12.61 For any a £ (0, 1), let us observe that, for any s > 0, 

(5 ' 4) ^~ S) = (l + snas + l) 

(5.5) = / (1 - e- s y)u a (y)dy, 

Jo 

where u a {y) = r( - 1 _ ct e ^ 1 ^ e -( ) / Q )a+i • We easily check that u a (y)dy € V and hence ^> is a Laplace 
exponent of a Levy process which drifts to —00. Next, we know, see e.g. [25], that 

T j Q — Ot 

H~ ~ a ' 
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where H is the negative of the subordinator having Laplace exponent 

~ , . aT(as + 1) 
T(a(s - 1) + 1) 

Observing that <p—(— s) = s + 1), we deduce, from (15. 2ft . that 

(5.6) ?n H - (x) = ^^<fe (a~ 1/a ) , x > 0, 

from which we readily get the expression t|2.10|) . Then, we recall the following power series 
representation of positive stable laws, see e.g. [291 Formula (14.31)], 



oo 



9. 



(-1)" 

'—an)n\ 



-(l+an) X>Q _ 



n=l 

a-1. 



Then, by means of an argument of dominated convergence justified by the condition lim s _ 
0, we get, for all x > 0, that 

= fe+ nfc=i^+(- fc ) ;c n i 
a ^— ' r(— an)n! 

n=l 

where we used the identities (|4.8p . E[— Y"i] = /c + and ip+(— = ~~ ^+( — k). The fact that the 
series is absolutely convergent is justified by using classical criteria combined with the Euler's 
reflection formula T(l — z)T(z) sin(-7rz) = ir with the asymptotics 

(5.7) Y^^j = ^ (1 + ° aS Z ^ °°' |ar5(z)! < 7F ' 

see e.g. [191 Chap. 1]. We complete the proof by mentioning that Simon [30] proved recently 
that the positive stable laws are MSU if and only if a < 1/2 which implies, from (|5.6p . that 
is also MSU in this case. 
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